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NOMENCLATURE 


A . Area  of  bearing  surface,  in^ 

A, B... Linear  regression  coefficients 

b . Width  of  flow,  in 

C . Fluid  specific  heat,  Btu/lbs*°F 

C^. . . .Constant  of  integration 
Cj ....  Constant  of  integration 
Ef ....  Frictional  energy  loss,  in*lbs/sec 
Ep....Pvimp  energy  loss,  in*lbs/sec 

F . Force,  lbs 

h . Film  thickness,  in 

J . Conversion  factor,  ft*lbs/Btu 

K . Power  Law  constant,  Pa*sec" 

Kg.  .  .  .Bingham  and  Hyperbolic  model  constant.  Parsec 

M . Frictional  torque,  in*lbs 

N . Angular  velocity,  RPM 

n . Power  Law  exponent 

P . Pressure,  Ibs/in^ 

P . Supply  pressure,  Ibs/in* 

AP. ...  Pressure  difference  across  bearing,  Ibs/in^ 

Q . Flow  rate,  in^/sec 

R . Outside  radius,  in 

Ro....Step  radius,  in 

Rj. . . .Thermister  one  resistance,  kOhms 

Rj. . . .Thermister  two  resistance,  kOhms 

v7i 

R . Reynolds  number,  — ,  Nondimensional 

r . Radius  value  between  R„  and  R,  in 

T . Temperature,  Rankine 

Tq. . . .Baseline  temperature,  Rankine 
T^. . . .Thermister  one  temperature,  Rankine 
Tj .... Thermister  two  temperature,  Rankine 
At. . . .Temperature  difference,  °F 

V  . Fluid  velocity,  in/sec 

v(y) . .Fluid  velocity,  in /sec 

9/j  ... 

....  Partial  derivative  of  film  thickness 
dyi 

Y  . Fluid  weight  density,  Ibs/in’ 
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Y  . Power  Law  model  shear  rate,  dv/dy 

p, . Fluid  viscosity,  reyns 

Oy. . . .Standard  deviation  along  y-axis 
T . Shear  stress,  Ibs/in^ 

Xg. . . .Bingham  and  Hyperbolic  model  yield  point.  Pa 
Xn-  .  .  .Newtonian  shear  stress,  Ibs/in^ 

■^Non-N-  .Non -Newtonian  shear  stress,  Ibs/in^ 

V  . Kinematic  fluid  viscosity,  in^/sec 

(0 . Angular  velocity,  rad/sec 

cOh.  .  .  .Uncertainty 


Abstract  of  Thesis  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of 
the  Requirements  for  the  Degree  of  Master  of  Science 

NON-NEWTONIAN  TEMPERATURE  AND  PRESSURE  EFFECTS 
OF  A  POWDER  LUBRICANT  SLURRY  IN  A 
ROTATING  HYDROSTATIC  STEP  BEARING 

By 

William  Edward  Finn 
May  1993 

Chairperson:  Dr.  J.E.  Peterson 

Major  Department:  Mechanical  Engineering 

The  purpose  of  this  research  was  to  investigate  the 
pressure  and  temperature  effects  of  graphite  powder 
lubricant  added  to  a  carrier  fluid  of  ethylene  glycol  as 
applied  in  a  rotating  hydrostatic  step  bearing.  After  the 
testing  apparatus  was  designed  and  modified,  the  temperature 
and  pressure  profiles  were  determined  analytically  and 
experimentally . 

The  rheological  behavior  of  the  non-Newtonian  lubricant 
was  modeled  using  a  Power  Law  model  developed  previously. 
This  model  has  been  proven  to  match  experimental  data  and 
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was  therefore  the  tool  used  for  comparison  with  research 
data. 

Ethylene  glycol  was  used  as  the  Newtonian  lubricant  for 
the  initial  experimentation,  which  provided  a  check  on  the 
test  apparatus  and  a  comparison  with  the  non-Newtonian 
graphite  slurry. 

The  measured  data  for  the  temperature  increase 
dependency  on  rotation  for  both  fluids  compared  favorably  to 
the  mathematical  predictions.  A  significantly  higher 
temperature  difference  was  seen  in  the  non-Newtonian 
lubricant  due  to  the  higher  shear  rates. 

The  pressure  profile,  while  not  directly  dependent  on 
rotation  as  shown  in  the  mathematical  model,  demonstrated  a 
reduction  based  upon  higher  rotation  speeds.  This  reduction 
was  greater  for  the  non-Newtonian  lubricant.  Due  to  the 
higher  temperatures  seen  in  the  non-Newtonian  lubricant  and 
temperature  dependent  Power  Law  constants,  a  pressure  loss 
is  expected.  For  the  Newtonian  lubricant,  viscosity  is 
dependent  upon  fluid  temperature. 

Therefore,  the  effects  of  operating  speed  and 
temperature  on  a  non-Newtonian  lubricant  should  be 
considered,  as  well  as  their  greater  load- carrying  capacity. 


viii 


CHAPTER  1 

BACKGROUND  AND  SIGNIFICANCE  OF  RESEARCH 

Background  of  Hydrostatic  Bearings 
Hydrostatic  bearings  were  first  demonstrated  in  I85i  by 
L.D.  Girard  through  an  invention  involving  high  pressure 
water- fed  bearings  used  for  a  railway  propulsion  system  and 
patented  in  1865  (i,  p.  225;  2,  p.  435) .  This  externally 
pressurized  or  hydrostatic  system  was  further  demonstrated 
at  the  1878  Paris  Industrial  Exposition  when  a  heavy  block 
with  four  feet  was  floated  on  a  thin  film  of  lubricant  that 
was  pumped  down  each  leg  (i,  2) .  The  first  analysis  of  a 
thrust  bearing  was  provided  by  Lord  Rayleigh  in  1917  using 
the  hydrostatic  system  and  solving  the  load,  flowrate  and 
frictional  torque  (1;  2,  p.  436) . 

Current  applications  include  the  analysing  magnet  of 
the  Nuclear  Structure  Facility  at  the  SERC  Daresbury 
Laboratory  (i,  p.  227) ,  the  Halle  optical  telescope  (2,  p. 
436) ,  the  2 10 -foot -diameter  tracking  antenna  at  the  NASA-Jet 
Propulsion  Laboratory  Deep  Space  Instrumentation  Facility 
(3,  p.  33-51),  and  Denver’s  Mile  High  Stadium  (4,  p.  78) . 
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All  applications  involve  moving  extreme  weight  with  minimal 
effort  and  in  some  cases,  high  precision.  In  1972,  the 
Institution  of  Mechanical  Engineers  published  proceedings  on 
a  conference  solely  devoted  to  externally  pressurized 
bearings  (2,  p.  437) .  Since  that  time,  more  effort  has  been 
spent  by  industry  on  the  study  of  the  lubricants  as  most  of 
the  basic  mechanical  principles  of  the  hydrostatic  bearing 
had  been  solved. 

Significance  of  the  Research 
The  addition  of  solid  lubricants  to  Newtonian  fluids 
has  been  the  subject  of  extensive  research.  One  major 
reason  for  the  addition  of  solid  lubricants  is  the  greater 
stability  of  these  mixtures  when  subjected  to  extreme 
temperatures.  The  aerospace  industry,  military,  and  NASA 
all  report  using  solid  lubricants  because  of  their  greater 
load  and  temperature  capabilities  (5,  p.  5;  6,  p.  5,  83-86) . 
Many  comme : cial  applications  exist  for  a  Graphite/Glycol 
lubricant  mixture  including  switch  gear  equipment  and  high 
temperature  oven  bearings  (5,  p.  85-86). 

The  U.S.  gas  turbine  industry  has  started  a  new  program 
in  cooperation  with  NASA  and  the  Department  of  Defense 
called  the  Integrated  High  Performance  Turbine  Engine 
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Technology  (IHPTET)  program  (7,  p.  114) .  The  desired  double 
propulsion  capability  from  the  large  speeds  in  the  gas 
turbine  shafts  result  in  high  temperatures  which  have 
further  emphasized  the  requirement  for  nontraditional 
lubricants.  The  approximate  working  temperatures  of 
different  types  of  lubricants  are  as  follows  (8,  p.  2) : 
Petroleum  Lubricants  below  100°C 

Synthetic  Lubricar  v,s  200  to  315°C 

Powder  Lubricants  425  to  675°C 

The  powder  lubricant  is  supplied  to  the  system  through 
pressurized  air  or  slurries  containing  the  powder  (7,  p. 

114) .  These  non- Newtonian  lubricants  are  the  subject  of 
current  research  towards  more  accurate  behavior  prediction. 


CHAPTER  2 

OBJECTIVE  OF  THE  RESEARCH 


This  research  investigates  the  pressure  and  temperature 
affects  of  graphite  powder  lubricant  added  to  a  carrier 
fluid  of  ethylene  glycol  as  applied  in  a  rotating 
hydrostatic  step  bearing.  The  graphite  powder  was  provided 
by  Union  Carbide  and  the  ethylene  glycol  was  provided  by 
Fisher  Scientific.  The  methodology  included  laboratoij' 
testing  and  mathematical  modeling  for  comparison  with  the 
experimental  results. 

An  existing  laboratory  test  rig  was  extensively 
modified  to  impart  an  angular  velocity  to  a  hydrostatic  step 
bearing  and  to  measure  temperature,  pressure,  flowrate, 
revolutions  per  minute  (RPM)  ,  and  film  thickness. 
Measurements  for  the  ethylene  glycol  were  obtained  for  use 
as  a  check  on  the  testing  rig  and  for  establishing  baseline 
data. 

The  goal  was  to  formulate  temperature  equations 
dependent  upon  the  rotational  velocity  and  compare  this 
affect  with  the  pressure  distribution  across  the  bearing 
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radius.  Potential  inertial  affects  from  the  step  bearing's 
rotation  and  its  effects  on  the  graphite  particles  in  the 
ethylene  glycol  would  also  be  examined  in  the  pressure 
profile . 

An  understanding  of  these  fundamental  lubrication 
principles  in  a  non-Newtonian  fluid  is  necessary  for 
developing  a  powder  lubricant  slurry  approach  in  future 
applications . 


CHAPTER  3 

RESEARCH  EQUIPMENT  AND  METHOD 

Laboratory  Test  Equipment 

The  laboratory  testing  was  an  important  portion  of  this 
research.  Significant  time  and  effort  was  spent  designing, 
fabricating,  and  redesigning  the  modifications  to  the 
existing  test  rig.  The  modified  hydrostatic  test  rig 
(Figure  3-1)  contains  a  flat  base  and  a  parallel  hydrostatic 
step  bearing  (items  5  and  6,  Figure  3-1) .  The  bearing 
consisted  of  an  upper  surface  with  an  outside  diameter  of 
four  inches  and  an  inside  diameter  of  two  inches.  It  was 
connected  to  the  shaft  and  gear  assembly  through  a  ball 
bearing  assembly  designed  and  fabricated  to  maintain  a 
parallel  surface  to  the  base  while  rotating.  The  resultant 
constant  film  thickness  was  measured  by  three  LVDT 
displacement  sensors  (item  4,  Figure  3-1)  mounted  on 
supports  which  were  placed  120  degrees  apart  on  top  of  the 
hydrostatic  step  bearing.  Since  three  points  define  a 
plane,  the  use  of  three  sensors  was  sufficient.  The  voltage 


6 


8 


reading  of  the  LVDT  displacement  sensors  was  compared  to  the 
calibration  and  provided  film  thickness.  The  specifications 

for  the  LVDT  sensors  are  as  follows: 

Make  Trans -Tek,  inc. 

Model  DC-DC  Series  240 

02141-0000 

Input  Volts  DC  (v)  6  to  30 

Linearity  %  .  5 

Freq.  Response  3dB  down  140  Hz 

Temperature  Range  -50  C  to  120  C 

Resolution  Infinite 

The  film  thickness  could  be  adjusted  through  the  fine 
adjusting  thread  (item  l.  Figure  3-1)  that  connected  the 
screw  shaft  to  the  frame  (item  2,  Figure  3-1) .  The  gear 
assembly  was  belt  driven  by  a  controllable,  variable  speed 
motor  (item  3,  Figure  3-1) .  The  RPM  was  measured  directly 
through  a  Pioneer  Digital  Photo  tachometer,  model  DT-36. 


The  motor  specifications  are  as  follows: 


Make 

Model 

HP 

RPM 

Armature  Volts 
Armature  Current 


Dayton  Permanent 
Magnet  DC 
2M168B 
1/2 
1725 
90  VDC 
5 . 5  amps 


Two  thermistors  were  installed  in  the  base  of  the 


bearing  surface  (item  7,  Figure  3-1)  and  measured  the 
temperature  of  the  lubricant  entering  the  recess  pocket  and 
at  the  perimeter  of  the  hydrostatic  step  bearing.  The 
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measured  resistance  was  compared  to  the  thermister's 
calibration  and  provided  the  temperature. 

The  lubricant  was  supplied  to  the  bearing  surfaces  by 
an  external  pump  through  an  internal  orifice.  The  pump's 
flow  rate  was  adjustable.  It's  specifications  are  as 
follows : 

Make  Micropump  Corp. 

Model  415 

Motor/Pump  Coupling  Magnetic  Drive 

Maximum  Flow  Rate  630  ml /min 

Rated  Maximum  Pressure  20  psi 

The  fluid  flow  was  routed  through  a  rotameter  flow 
meter.  Comparing  the  reading  with  the  flow  meter's 
calibration  provided  the  flow  rate. 

The  fluid  film  pressures  were  measured  with  dial 
pressure  gages  at  four  locations  across  the  bearing  surface 
located: 

1)  within  the  recess  pocket 

2)  1/4  inch  from  the  inside  pocket 

3)  1/2  inch  from  the  inside  pocket 

4)  3/4  inch  from  the  inside  pocket 

Research  Method 

The  experimental  test  rig  was  used  to  determine  the 
effects  of  powder  lubricants  in  a  carrier  fluid  subjected  to 
rotation  in  a  hydrostatic  step  bearing.  The  testing  method 
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is  described  in  Chapter  5.  First,  the  various  measuring 
devices  were  calibrated. 

Uncertainty  Analysis  and  Calibration 
The  uncertainty  analysis  for  the  flowmeter,  thermister, 
and  LVDT  data  was  calculated  using  a  least  squares  fit 
technique  to  determine  the  deviation  width.  Since  each 
cuirve  fit  was  linear,  the  formula  (9,  p.  158)  used  to 
determine  the  standard  deviation  is  written  as  equation 
(3-1)  where  the  coefficients  A  and  B  are  determined  through 
the  linear  regression  curve  fit.  Further,  was  multiplied 
by  1.7  in  order  to  obtain  10:1  odds  or  a  90%  probability 
that  a  measurement  will  fall  within  1.7  standard  deviations. 

Cij  =  j^'Liyi-A-Bxif  (3-1) 

The  rotameter  flowmeter  was  calibrated  by  adjusting  the 
flowrate  through  the  pump,  marking  the  flowmeter  reading, 
and  timing  the  fill  time  of  the  liquid  into  a  beaker  of 
known  volume.  A  linear  regression  or  least-squares  fit  (9, 
p.  154)  for  a  line  of  the  measured  data  provided  the 


calibration  equation: 
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y  =  5. 55jc- 205.69  (3-2) 

where  y  is  the  flow  rate  in  ml/min  and  x  is  the 
flowmeter  reading.  The  calibration  curve  is  shown  in  figure 
3-2.  The  flowmeter  uncertainty  calculated  through  equation 
(3-1)  for  both  liquids  was  ±11.8  ml/min. 


*  Data  - Eq.  (3-2) 

Figure  3-2.  Flowmeter  calibration  curve 

The  thermisters  were  calibrated  by  placing  them  in  a 
water  bath,  adjusting  the  temperature,  reading  the  mercury 
thermometer  temperature,  and  recording  the  thermisters' 
resistance.  The  uncertainty  of  the  thermometer  using  half 
the  least  count  was  ±0.1°F.  Thermisters  function  on  a 
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principle  that  their  resistance  is  lower  at  higher 
temperatures.  This  dependency  is  logarithmic  and  shown  by 
the  equation: 


where  and  are  the  baseline  resistance  in  kOhms 
and  temperature  in  Rankine,  respectively.  The  constants  A 
and  B  were  determined  through  calibration  for  the  two 
thermisters  using  linear  regression.  Rearranging  the 
results  and  solving  for  resistance  as  a  function  of 
temperature  provided  the  calibration  equations: 


«2 


(3-4) 

(3-5) 


where  R^,  T^  and  Rj,  T^  are  the  resistances  and 
temperatures  of  thermisters  one  and  two,  respectively.  The 
calibration  curves  are  shown  in  figures  3-3  and  3-4.  Their 
uncertainty  calculated  through  equation  (3-1)  was  ±  l  Ohm 


(.001  kOhms )  . 


The  LVDTs  were  calibrated  by  unscrewing  the  shaft  of 
known  adjusting  thread  through  equidistant  intervals 
beginning  at  zero  when  the  hydrostatic  step  bearing  was 
touching  the  flat  base.  The  voltage  reading  of  each  LVDT 
was  recorded  at  these  intervals.  Linear  regression  of  the 
measured  data  for  each  LVDT  provided  the  following 
calibration  equations: 

LVDT  A:  = -1 48.4jc^  +  1 9.44  (3-6) 

LVDT  B:  = -133.8X5  +  19.78  (3-7) 

LVDT  c:  yc  = -97. 6xc+ 19.97  (3-8) 

where  y  is  the  LVDT  reading  in  volts  and  x  is  the  film 

thickness  in  inches.  The  calibration  curves  are  shown  in 
figure  3-5.  The  uncertainty  calculated  through  equation 
(3-1)  was  ±0.04  volts  for  LVDT  A,  ±0.07  volts  for  LVDT  B, 
and  ±0.05  volts  for  LVDT  C. 

The  uncertainty  for  the  tachometer  was  ±i  RPM  as  given 
by  the  manufacturer.  It  was  verified  by  repeatedly  setting 
the  maximum  RPM  value  and  noting  that  the  value  was 
consistently  within  ±l  of  670  RPM.  The  uncertainty  for  the 
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Figure  3-5.  LVDT  calibration  curves 


dial  pressure  gages  was  ±0.5  psi  using  half  the  least  count. 

The  film  thickness,  h,  was  calculated  as  an  average  of 
the  three  LVDT  readings,  in  order  to  determine  the  overall 
uncertainty  of  the  film  thickness  during  the  testing, 
equations  (3-6),  (3-7),  and  (3-8)  were  rearranged  to  solve 

for  film  thickness  based  upon  LVDT  voltage.  The  three 
equations  were  averaged  resulting  in 


h  =  =  -.0022yA  -  .0025yB- .0033yc  + .161  ( 3 - 9 1 


16 


A  root  -  sum- square  method  (9,  p.  73;  10,  p.  15;  11,  p. 

254)  was  employed  as  shown  in  equation  (3-10)  where  h  is 

3/2 

determined  from  equation  (3-9)  and  ^  is  its  partial 
derivative  with  respect  to  LYDT  A,  B,  or  C.  The  uncertainty 
of  LVDT  A,  B,  or  C  is  co^  as  derived  earlier. 


0)/,= 


(3-10) 


The  overall  uncertainty  of  h,  the  film  thickness,  was 
±0.0003  inches.  This  is  important  because  the  film 
thickness  in  the  mathematical  modeling  of  the  pressure 
profile  is  raised  to  the  third  power  resulting  in  a  three 
times  greater  uncertainty. 


CHAPTER  4 

MATHEMATICAL  MODELING 


The  mathematical  calculations  this  chapter  were 
formulated  in  order  to  model  the  temperature  change  and  the 
pressure  profile  in  a  hydrostatic  step  bearing  while  varying 
the  film  thickness  and  revolutions  per  minute  (RPM) .  The 
comparisons  of  the  measured  data  with  the  mathematical 
modeling  are  discussed  in  Chapter  5. 


Figure  4-1.  Schematic  of  Hydrostatic  Step  Bearing 
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The  equations  derived  for  the  temperature  and  pressure 
relationships  of  the  hydrostatic  step  bearing  shown  in 
figure  4-l  involve  some  basic  assumptions  (4,  pp.  80-88;  8, 
p.  10)  ; 

1.  The  flow  is  laminar.  The  maximum  Reynolds  number 
using  the  greatest  flow  rate  and  the  radius  at  the 
step  is  25.9.  A  Reynolds  number  less  than  1000  in 
bearing  applications  is  assumed  to  be  laminar  flow 
(4,  p,  63;  12,  p.  343-346) .  The  flow  is  pressure 
induced. 

2 .  The  inertial  forces  resulting  from  the  rotation  of 
the  bearing  are  negligible  when  compared  to  the 
viscous  shear  forces  and  can  be  neglected. 

3.  The  fluids  are  incompressible. 

4.  The  pressure  is  a  function  of  the  radius, 
neglecting  any  side  leakage. 

5.  The  temperature  is  constant  across  the  film 
thickness . 

6.  All  the  heat  generated  appears  in  the  fluid  and 
corresponds  to  its  temperature  rise. 

7.  The  pressure  at  the  outside  of  the  bearing  is 
assumed  to  be  ambient. 
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Temperature  Equations 

Because  of  the  different  properties  of  Newtonian  and 
non-Newtonian  fluids,  two  different  mathematical  model 
relations  were  used. 

Ethylene  Glycol  Relations 

Pure  ethylene  glycol  is  a  Newtonian  fluid  and 
therefore,  the  following  relations  apply  according  to  the 
geometry  of  the  test  rig  as  shown  in  the  schematic  (figure 
4-1)  . 

The  pump  loss  is  given  by  the  relation  (4,  p.  86) , 

Ep  =  QAP  (4-1) 

where  Q  is  the  flow  rate  in  cubic  inches  per  second, 

AP  is  the  pressure  difference  in  psi,  and  Ep  is  energy  in 
inch-pounds  per  second.  The  friction  loss  on  the  bearing 
due  to  rotation  can  be  solved  beginning  with  Newton's  fluid 
friction  law. 


(4-2) 
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Referring  to  Figure  4-1  and  substituting  the  film 
thickness,  h,  for  the  differential  width,  this  is  equivalent 
to  Newton's  equation  of  viscous  friction  drag. 


F  =  n/ll 


(4-3) 


Substituting  jir^  for  A  and  considering  a  differential 
annulus  results. 


dF=\ll{2Tzrdr) 


(4-4) 


When  one  of  the  contacting  surfaces  is  stationary, 
v=0)r,  where  co  is  the  angular  velocity  in  radians  per 
second.  Multiplying  both  sides  by  r  to  determine  the  torque 
on  an  elemental  annulus  of  radius  r  and  width  dr  is 


rdF  =  dM  =  ^^^r^dr 

h 


(4-5) 


The  total  frictional  torque  (4,  p.  88)  becomes 


2n[i(a  rR 


iRo  h  \_  A  _ 


(4-6) 
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where  n  is  the  viscosity  in  reyns,  co  is  the  angular 
velocity  in  radians  per  second,  h  is  the  film  thickness  in 
inches,  is  the  radius  to  the  step  in  inches,  and  R  is  the 
radius  to  the  outside  edge  of  the  bearing  in  inches. 

The  friction  loss  is  given  by  the  equation  (4,  p.  88) 

(4-7) 

where  is  energy  in  inch-pounds  per  second,  N  is  the 
rotation  in  revolutions  per  minute  (RPM) ,  and  M  is  the 
torque  in  inch-pounds.  Therefore,  the  total  energy  loss 
across  the  bearing  from  r=o  to  r=R  can  be  expressed  as 
follows  (4,  p.  90)  , 

Ej-  =  Ep+Ef=\2JQyCAT  (4-8) 

where  12 J  is  the  conversion  factor  for  foot-pounds  to 
British  thermal  units  (Btu) ,  Q  is  the  flow  rate  of  the  fluid 
in  cubic  inches  per  second,  y  is  the  weight  density  of  the 
fluid  in  pounds  per  cubic  inches,  C  is  the  specific  heat  of 
the  fluid  in  Btu  per  pounds- °F,  and  AT  is  the  temperature 
difference  across  the  bearing  in  °F.  The  temperature 


Ef  = 


2nMN 

60 
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difference  can  be  solved  for  as  a  function  of  the  rotation 
since  the  flow  rate,  pressure  difference,  and  film  thickness 
are  held  constant.  The  pressure  difference  from  r=0  to  r=R 
equals  the  supply  pressure  since  the  pressure  at  the  outside 
of  the  bearing,  P(R) ,  is  assumed  to  be  ambient.  Therefore, 
substituting  the  values  for  and  Ep  and  rearranging  to 
solve  for  AT  results 


Ar= 


limOJQiCh 


(4-9) 


Graphite /Ethylene  Glycol  Relations 

A  mixture  of  ethylene  glycol  and  graphite  powder  is  a 
non-Newtonian  fluid  and  requires  different  fundamental 
modeling.  Currently,  there  are  three  accepted  rheological 
models  (7,  pp.  115-116;  8,  p.  15) .  First,  the  Bingham  model 
is  stated: 

'^Non-N  —  '^0+  Koy  ( 4  - 10 ) 

where  z^on-ti  is  the  non-Newtonian  shear  stress  in 
Pascals,  To  is  the  yield  point  in  Pascals,  Ko  is  the  Bingham 
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constant  in  Pascal -seconds,  and  y  is  the  shear  rate  in 
seconds'^.  This  model  has  shown  to  deviate  significantly 
from  measurements  at  low  shear  rates,  however,  it  is 
accurate  at  high  shear  rates  where  many  lubricant  films 
operate  (7,  p.  115) .  Second,  the  Hyperbolic  model  is 
stated: 


'lNon-N  =  [{Ko'if  +2Ko'lo'{^’^  (4-11) 

where  the  shear  stress,  yield  point,  Bingham  constant, 
and  shear  rate  variables  are  the  same  as  the  Bingham  model. 
This  function  is  similar  to  the  Bingham  model  at  high  shear 
rates  and  has  proven  to  match  better  with  results  at  low 
shear  rates  (7,  p.  116) . 

Finally,  the  Power  Law  model  is  stated: 

Non-N  ~  (4-12) 

where  fNon-N  is  the  non-Newtonian  shear  stress  in 
Pascals,  K  is  the  Power  Law  constant  in  Pascal -seconds",  y 
is  the  shear  rate  in  seconds",  and  n  is  a  nondimens ional 
Power  Law  constant.  This  model  has  proven  useful  for 


describing  shear  stress -shear  rate  relations  for 


non-Newtonian  fluids  without  a  yield  point  and  is  useful  for 
fluid  flow  formulations  (7,  p.  116) . 

Since  it  has  been  shown  that  the  Power  Law  model  can 
reasonably  model  a  1:8  ratio  of  graphite  powder  to  ethylene 
glycol  (7,  p.  116;  8,  pp.  31-32;  13,  pp.  86-87)  and  the 
mathematical  expression  is  uncomplicated  to  evaluate,  it  was 
used  in  this  research. 


Table  4-1.  Power  Law  Constants 


Graphite -Ethylene  Glycol  (1:8) 

Temp,  °C 

K,  Pa- sec" 

n 

20 

0.0572 

0.973 

30 

0.0448 

0.971 

40 

0.039 

0.969 

50 

0.0361 

0.969 

The  Power  Law  constants  used  were  taken  from  Table  4-1 
(7,  p.  116)  dependent  on  the  temperature  of  the  fluid. 

The  pump  loss  given  by  equation  (4-1)  is  still 


applicable.  The  friction  loss  on  the  bearing  due  to 
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rotation  can  be  determined  beginning  with  the  Power  Law 
model,  equation  (4-12),  rewritten  as  follows: 

X=k[^J  (4-13, 

This  is  equivalent  to  the  equation  of  viscous  friction 

drag. 


F  =  KA  I 

h 


(4-14) 


Substituting  for  A  and  considering  a  differential 


annulus  results. 


dF  =  K\jX  (Inrdr) 


(4-15) 


Since  one  of  the  contacting  surfaces  is  stationary, 
v=0)r,  where  O)  is  the  angular  velocity  in  radians  per 
second.  Multiplying  both  sides  by  r  to  determine  the  torque 
on  an  elemental  annulus  of  radius  r  and  width  dr  results. 


rdF=dM=  2kk{^X  (r^^^dr) 


(4-16) 


The  total  frictional  torque  becomes 


M=2%K 


This  is  equivalent  to  equation  (4-6)  for  the  Newtonian 
fluid  where  n=l  and  K  is  the  viscosity.  A  conversion  factor 
is  necessary  since  the  Power  Law  constant,  K,  has  units  of 
Pascals  versus  pounds  per  square  inch  (psi) .  The  friction 
loss  from  equation  (4-7)  expands  to  the  following: 


p  -  n^KNf  Vr 

15  [30h  J  I  n+3 


(4-18) 


Therefore,  the  total  energy  loss  from  r=0  to  r=R  can  be 
determined  according  to  equation  (4-8) .  The  temperature 
difference  is  solved  as  a  function  of  the  rotation  since  the 
flow  rate,  pressure  difference,  and  film  thickness  are  held 
constant.  The  pressure  difference  from  r=0  to  r=R  equals 
the  supply  pressure  since  the  pressure  at  the  outside  of  the 


bearing,  P(R) ,  is  assumed  to  be  ambient.  Therefore, 


derivations  of  equations  (4-2)  and  (4-12) .  The  fluid 
flow  in  the  hydrostatic  step  bearing  is  pressure  induced. 
Ethylene  Glycol  Relations 

For  the  geometry  of  the  hydrostatic  step  bearing 
(Figure  4-1)  ,  the  pressure  distribution  on  a  unit  vol\ame  is 
shown  in  Figure  4-2. 

Summing  the  forces  in  the  horizontal  direction  and 
rearranging  is 


dx  dy 


(4-20) 


Substituting  equation  (4-2) ,  it  is  a  simplified  form  of 
the  Navier-Stokes  equations  (12,  p.343)  for  pressure  induced 
flow  between  two  parallel  surfaces, 


^  =  u— 

dx 


(4-21) 


Integrating  twice  across  the  y  direction,  assuming 


dp/dx  is  a  constant  across  the  film  thickness,  yields 
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Evaluating  the  constants  of  integration,  and  Cj, 
with  the  boundary  conditions  that  v(0)=0  and  v{h)=0  results 


''M  =  ->'*)  (4-23) 

Solving  for  the  flow  rate,  Q,  is  shown 

Q  =  \^Q^iy)dyb  (4-24) 

Integrating,  it  is  stated  (4,  p.  81;  7,  p.  118;  8,  p. 


^  \2\idx 


(4-25) 


where  Q  is  the  flow  rate  in  cubic  inches  per  second,  b 
is  the  plate  width  in  inches,  h  is  the  film  thickness  in 
inches,  ji  is  the  viscosity  in  reyns,  and  dp/dx  is  the 
pressure  gradient.  Equation  (4-25)  is  commonly  used  in 
radial  flow  situations.  Since  ethylene  glycol  is  a 
Newtonian  fluid.  Replacing  b  by  2Kr,  dx  by  dr,  integrating 
and  solving  for  the  pressure  at  any  radius,  r,  becomes  (4, 
p.  82;  8,  p.l4;  14,  p.  7) 
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I 

I 

I 


P  = 


izh^ 


(4-26) 


where  P  is  the  pressure  at  any  point  along  the  radius 
in  pounds  per  square  inch  (psi) ,  R  is  the  outside  radius  of 
the  bearing,  r  is  any  radius  between  R„  and  R,  Q  is  the  flow 
rate  in  cubic  inches  per  second,  h  is  the  film  thickness  in 
inches,  and  p,  is  the  viscosity  in  reyns. 

Graphite /Ethylene  Glycol  Relations 

Using  the  same  assumptions  as  the  ethylene  glycol 
relations  with  the  exception  of  equation  (4-12)  for  the 
Power  Law  Model  described  earlier,  it  is  stated  (7,  p.  118; 
8,  p.  16;  14,  p.  7) 


2rH-\ 

n 


^2«+l ^ 


(4-27) 


Replacing  b  by  27rr,  dx  by  dr,  integrating  and  solving 
for  the  pressure  at  any  radius,  r,  becomes  (8,  p.  16;  14,  p. 
8) 


/  \n  /  \ 


1-n 


A 

J 


2 

h 


2m+1 


(4-28) 
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The  Power  Law  constants,  K  and  n,  are  determined  from 
table  4-1  for  equations  (4-22)  and  (4-23)  dependent  upon  the 
temperature  of  the  fluid. 

For  this  research,  the  pressure  /alues  for  both  fluids 
were  nondimensionalized  by  dividing  hy  the  inlet  supply 
pressure.  The  flow  rate  and  film  thickness  were  constant  so 
that  the  RPM  could  be  adjusted  in  order  to  observe  the 
rotational  effects  on  the  pressure  at  the  set  points  along 
the  radius. 


CHAPTER  5 

EXPERIMENTAL  RESULTS 


Experimencal  Method 

The  initial  testing  with  the  laboratory  test  rig  was 
performed  with  pure  ethylene  glycol  as  the  lubricant.  This 
data  provided  a  check  on  the  instruments  and  equipment  and  a 
comparison  for  studying  the  effects  of  adding  graphite 
powder  to  the  ethylene  glycol. 

The  film  thickness  was  set  based  upon  the  LVDT 
readings.  The  flow  rate  of  the  lubricant  was  adjusted  so 
that  the  maximxim  pressure  of  the  dial  gages  would  not  be 
exceeded.  The  motor  was  energized  and  readings  were  taken  of 
the  LVDT  voltages,  flowmeter,  pressure  gages,  and 
thermisters  as  the  rotation  was  adjusted  from  lOO  to  670 
RPM.  When  the  Graphite/Ethylene  Glycol  slurry  was  the 
lubricant,  equivalent  film  thicknesses  as  the  pure  ethylene 
glycol  were  obtained  and  the  above  procedure  was  performed. 

Temperature  Results 

The  purpose  of  the  thermister  data  was  to  detenr.xne  the 
relation  between  the  temperature  difference  at  the  inlet  and 
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the  edge  of  the  bearing  and  the  rotation  speed  (RPM)  .  This 
experimental  data  were  then  compared  with  the  mathematical 
predictions  as  developed  in  Chapter  4 .  The  data  shown  in 
Figures  5-1  and  5-2  closely  match  the  predictions  within  the 
uncertainty.  The  repeatability  was  within  experimental 
uncertainty. 
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D  h=  0,002 

»  h= 0.003 

♦  h  =  0.004 

Eq.  (4-9) 

.  Eq.  (4-9) 

— -  Eq,  (4-9) 

Figure  5-1.  Ethylene  Glycol  Temperature  Profile 
Flow  rate,  h=0.002:  77  ml/min 

h=0.003:  153  ml/min 

h=0.004;  371  ml/min 

Supply  Pressure,  h=0.002:  40  psi 

h=0.003:  25  psi 

h=0.004;  25  psi 


RPM 


Figure  5-2  clearly  shows  that  the  temperature 
difference  at  higher  RPM  for  the  graphite  slurry  is  almost 
double  the  temperature  difference  shown  in  Figure  5-1  at 


^  ...#•■ 

T' 


D  h=.002  *  h=.003  X  h=.0O4 

-  Eq.  (4-19)  .  Eq.  (4-19)  . Eq.  (4-19) 


Figure  5-2. 


Graphite/Ethylene  Glycol  Temperature 
Profile 

♦  Flow  rate,  h=0.002:  46  ml/min 

h=0.003:  114  ml/min 

h=0.004;  136  ml/min 

♦  Supply  Pressure,  h=0.002:  40  psi 

h=0.003:  30  psi 

h=0.004;  15  psi 


equal  RPM  and  film  thickness  for  the  Newtonian  fluid, 
ethylene  glycol.  It  is  especially  prevalent  for  the  0.002 
inch  film  thickness.  This  demonstrates  that  the  presence  of 
the  graphite  powder,  thus  making  the  fluid  non -Newtonian, 
greatly  increases  the  temperature  change.  For  a  rotating 


bearing  with  microstructures  in  the  lubricant,  the  increase 
in  temperature  and  subsequent  power  loss  is  expected  due  to 
an  increase  in  the  shear  stress  (15,  p.  106) .  As  discussed 
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in  Chapter  4,  the  Power  Law  constants  are  dependent  upon 
temperature.  The  values  used  for  the  rotation  speed  up  to 
500  RPM  corresponded  to  temperature  readings  within  the  30°C 
range.  When  the  rotation  exceeded  500  RPM,  the  temperature 
readings  were  higher  and  fell  within  the  40 °C  range. 
Therefore,  the  Power  Law  constants  were  adjusted  resulting 
in  the  discontinuity  shown  in  the  theoretical  curve  in 
Figure  5-2  at  h=0.002  inches.  The  uncertainty  for  the  data 
points  was  too  small  to  plot.  The  uncertainty  for  the 
mc.thematical  model  was  determined  as  discussed  in  Chapter  3 
with  equation  (3-2)  . 

Pressure  Results 

The  purpose  of  the  pressure  data  was  to  determine  the 
relationship  between  the  pressure  profile  and  the  rotation 
speed  (RPM)  .  The  inertial  and  temperature  effects  on  the 
pressure  are  assumed  to  be  negligible  when  compared  to  the 
shear  stress.  It  has  already  been  shown  that  the  addition 
of  the  graphite  particles  increases  the  pressure  at  lower 
flow  rates  thus,  increasing  the  load  capacity  of  the 
hydrostatic  step  bearing  (8,  p.  33;  14,  p.  16;  15,  p.  106; 
16,  pp.  32-33)  .  The  effect  of  rotation  on  the  pressure 
ratio  versus  th»-  radius  ratio  is  shown  in  Figures  5-3 
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through  5-8.  The  measured  data  compared  closely  with  the 
mathematical  predictions  as  developed  in  Chapter  4  at  the 
lower  rotation  speeds  within  the  uncertainty. 


Figure  5-3.  Ethylene  Glycol  Pressure  Profile 
♦  h=0.002;  Flow  rate=77  ml/min 


The  repeatability  was  within  experimental  uncertainty. 
The  uncertainty  for  the  data  points  was  too  small  to  plot. 
The  uncertainty  for  the  mathematical  model  was  determined  as 
discussed  in  Chapter  3  with  equation  (3-2) . 

The  inlet  pressure  varied  slightly  during  the  course  of 
the  experimentation  and  therefore,  the  pressures  were 


nondimensionalized  by  the  appropriate  supply  pressure 
reading  resulting  in  a  pressure  ratio  of  1  at  the  step  of 


Figure  5-4.  Ethylene  Glycol  Pressure  Profile 
♦  h=0.003;  Flow  rate=153  ml/min 


Figure  5-5.  Ethylene  Glycol  Pressure  Profile 


h=0.004;  Flow  rate=37i  ml/min 


0.9 

0.0 


temperature  of  the  bearing  increased  due  to  the  raised 
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frictional  loss  at  the  higher  rotations.  As  temperature 
increases,  the  viscosity  becomes  smaller  resulting  in  a 
smaller  pressure  ratio.  This  decrease  in  pressure  is  also 
indicated  in  Figures  5-4  and  5-5.  Preliminary  temperature 
measurements  of  the  interior  of  the  bearing  showed  higher 
temperatures  than  those  observed  at  the  inlet  or  outlet 
points  of  the  fluid  flow  emphasizing  the  presence  of  a 
greater  temperature  affect  on  the  interior  pressure  reading 
than  the  outlet  readings.  The  non-Newtonian  pressure 
behavior  shown  in  Figure  5-6  also  indicates  a  drop  in  the 
pressure  ratio.  This  decrease  in  the  pressure  ratio  is 
larger  than  shown  for  an  equal  film  thickness  and  RPM  of  the 
ethylene  glycol  in  Figure  5-3.  These  results  are  further 
supported  by  comparing  Figures  5-7  and  5-8  to  their  ethylene 
glycol  counterparts,  Figures  5-4  and  5-5.  As  shown  in  Table 

4- 1,  the  Power  Law  constants  for  the  graphite/ethylene 
glycol  slurry  are  also  temperature  dependent.  Therefore, 
the  non-Newtonian  graphite  slurry  is  also  temperature 
dependent.  The  graphite  slurry  demonstrated  a  higher 
pressure  drop  than  the  Newtonian  fluid,  ethylene  glycol, 
because  of  the  higher  temperatures  shown  in  Figures  5-i  and 

5- 2  and  discussed  previously. 
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The  decrease  in  pressure  is  significant  because  it 
translates  to  a  loss  of  load  carrying  capability  of  the 
bearing.  The  higher  load  carrying  capacity  of  the 
non-Newtonian  lubricant  is  one  of  the  reasons  for  their  use. 
The  pressure  is  still  higher  than  the  Newtonian  Ixibricant, 
but  it  must  be  considered  because  the  change  in  temperature 
experienced  at  extreme  rotations  will  be  substantially 
higher  than  Newtonian  lubricants  resulting  in  a  higher 
pressure  loss. 


CHAPTER  5 
CONCLUSIONS 


The  objective  of  this  research  was  to  determine  the 
effects  of  adding  a  powder  lubricant  to  a  Newtonian  carrier 
fluid  on  the  lubrication  temperature  and  pressure  profile. 
The  effects  were  shown  through  result  comparisons  between 
the  pure  ethylene  glycol  and  the  Graphite/Ethylene  Glycol 
slurry. 

A  detailed  uncertainty  analysis  of  the  measuring 
devices  was  undertaken  to  determine  deviations  between 
measured  and  predicted  data. 

Results  compared  favorably  with  the  mathematical 
models . 

The  data  presented  agreed  with  the  theoretical 
predictions  based  upon  the  rheological  properties  obtained 
from  previous  papers.  These  agreements  give  a  high  level  of 
confidence  for  the  use  of  the  equations  derived  in  this 
paper . 

It  was  noted  that  the  pressure  profile  of  the  graphite 
slurry  was  affected  more  by  the  higher  rotations  than  the 
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pure  Newtonian  fluid.  This  affect  can  be  expected  due  to 
the  higher  shear  rate  and  subsequent  temperature  increase. 
Since  viscosity  of  the  graphite  slurry  is  temperature 
dependent  and  the  potential  applications  involve  a 
significant  increase  to  the  rotation  speed,  it  is 
recommended  that  further  research  should  be  conducted  to 
study  the  interior  temperature  gains  and  their  affect  on  the 
pressure  profile  of  the  bearing.  Measuring  the  temperatures 
of  the  interior  of  the  bearing  in  conjunction  with  the  Power 
Law  constants'  curve  fit  as  functions  of  temperature  could 
result  in  a  complete  understanding  of  the  pressure  profile 
of  a  hydrostatic  step  bearing  subjected  to  a  rotation. 

Preliminary  testing  showed  that  the  temperature  inside 
the  bearing  is  higher  than  those  noted  at  the  inlet  or 
outlet  of  the  fluid  flow.  Therefore,  any  attempt  to 
mathematically  model  the  temperature  dependency  of  the 
viscosity  or  the  Power  Law  constants  would  have  been 
inaccurate . 

Another  factor  that  could  be  present  at  significantly 
higher  rotation  speeds  than  those  tested  in  this  research  is 
possible  inertial  effects  of  the  graphite  powder  in  the 
slurry  mixture.  That  effect  was  not  apparent  in  this 


research  since  an  increase  in  the  pressure  profile  would  be 
expected  for  such  an  occurrence. 

The  advantages  of  the  non-Newtonian  lubricant  are  their 
increased  viscosity  and  subsequent  higher  load  carrying 
capacity  and  their  ability  to  operate  at  higher 
temperatures.  However,  when  the  bearing  is  operated  at  high 
rotation  speeds,  the  higher  shear  rate  creates  an  increased 
temperature  which  could  have  a  weakening  affect  on  their 
pressure  profile.  This  affect  would  limit  their  load 
carrying  capacity,  although  it  still  exceeds  the  capacity  of 
the  Newtonian  lubricant.  This  thesis  has  shown  the 
fundamental  lubrication  principles  of  a  graphite  slurry 
non-Newtonian  fluid  that  was  pressure  induced  into  a 
rotational  hydrostatic  step  bearing.  Further  investigations 
into  higher  rotation  speeds  and  interior  temperature  gain 
are  recommended  to  more  fully  understand  the  non-Newtonian 
behavior  of  the  graphite/ethylene  glycol  lubricant. 
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APPENDIX  SAMPLE  UNCERTAINTY  CALCULATION 


Determining  the  uncertainty  in  the  film  thickness, 
equation  (3-9) ,  was  accomplished  through  the  use  of  equation 
(3-10)  which  is  expanded  as  follows: 


m  = 


2 

+ 


Substituting  the  values  of  the  partial  derivatives  of 
equation  (3-9)  and  the  uncertainties  of  each  LVDT  results: 

(0*  =  [(-.0022(±.04))^  +  (-.0025(±.07))^  +  (-.0033(±.05))^]  = 

03;,  =  ±0.00026  =  ±0.0003 

This  is  the  uncertainty  stated  on  page  16  and 
demonstrates  the  procedure  followed  for  all  the  uncertainty 
calculations . 
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